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Abstract — The problem of optimal control of power distri- 
bution systems is becoming increasingly compelling due to 
the progressive penetration of distributed energy resources in 
this specific layer of the electrical infrastructure. Distribution 
systems are, indeed, experiencing significant changes in terms of 
operation philosophies that are often based on optimal control 
strategies relying on the computation of linearized dependencies 
between controlled (e.g. voltages, frequency in case of islanding 
operation) and control variables (e.g. power injections, trans- 
formers tap positions). As the implementation of these strategies 
in real-time controllers imposes stringent time constraints, the 
derivation of analytical dependency between controlled and 
control variables becomes a non-trivial task to be solved. With 
reference to optimal voltage and power flow controls, this paper 
aims at providing an analytical derivation of node voltage and 
line current flows as a function of the nodal power injections 
and transformers tap-changers positions. Compared to other 
approaches presented in the literature, the one proposed here 
is based on the use of the [Y] compound matrix of a generic 
multi-phase radial unbalanced network. In order to estimate the 
computational benefits of the proposed approach, the relevant 
improvements are also quantified versus traditional methods. 
The validation of the proposed method is carried out by using 
both IEEE 13 and 34 node test feeders. The paper finally shows 
the use of the proposed method for the problem of optimal voltage 
control applied to the IEEE 34 node test feeder. 

Index Terms — Voltage/current sensitivity coefficients, unba- 
lanced electrical distribution networks, power systems optimal 
operation, smart grids. 



I. Introduction 

Optimal controls of power systems are often based on the 
solution of linear problems that link control variables 
to controlled quantities by means of sensitivity coefficients. 
Typical optimization problems refer to scheduling of gene- 
rators, voltage control, losses reduction, etc. So far, these 
categories of problems have been commonly investigated in 
the domain of high voltage transmission networks. However, 
during the past years, the increased penetration of distributed 
energy resources (DERs) in power distribution systems has 
raised the importance of developing optimal control strategies 
specifically applied to the operation of these networks (e.g. 1 1 1, 
0, 0, El, 0, [6 1). Within this context, it is worth noting 
that the solution of optimal problems becomes of interest only 
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if it meets the stringent time constraints required by real- 
time controls and imposed by the higher dynamics of these 
networks compared to the transmission ones. 

Typical examples of optimal controls that are not yet de- 
ployed in active distribution networks (ADNs) are voltage and 
power flow controls. Usually this category of problems has 
been addressed in the literature by means of linear-approaches 
applied to the dependency between voltages and power flows 
as a function of the power injections (e.g. iU, 0, 0, 0). 

The typical approach for the solution of this class of control 
problems is the use of the sensitivity coefficients through an 
updated Jacobian matrix derived from the load flow problem 
0, nni, HD, us, m. However, from the computational 
point of view, the main disadvantage of such a category of 
methods is that, for every change in the operation conditions 
of the network, an updated Jacobian matrix needs to be built on 
the basis of the network state and needs, then, to be inverted. 
This procedure involves non-trivial computation constraints for 
the implementation in real-time centralized or decentralized 
controllers. 

For this reason, the authors of [14J have proposed the 
direct computation of voltages and network losses sensitivity 
coefficients, based on the Gauss-Seidel formulation of the load 
flow problem, by making use of the [Z] matrix of a balanced 
network. Also, in |7| it has been proposed the use of the [Z] 
matrix along with the constant-current model for loads and 
generators. In |8] the sensitivity coefficients are proposed to 
be calculated starting from the network branch currents. The 
approach presented in (15), fT6\, fTf), fTSl, fT9l belongs to a 
class of methods typically derived from circuit theory and is 
based on the use of the so-called adjoint network. 

In order to increase the computational efficiency of this 
category of approaches, and to extend it to the inherent multi- 
phase unbalanced configuration of distribution networks, the 
main contribution of this paper is to provide a straightforward 
analytical derivation of node voltages and line currents sensi- 
tivities as a function of the power injections and transformers 
tap-changers positions. To this end, we propose to use the 
so-called [Y] compound matrix, which has the advantage of 
being sparse. 

Compared to [7| the approach here proposed takes into 
account the whole admittance matrix of the network. On 
the other hand the analytical derivation of sensitivities in 
was based on the approximated representation of the network 



lines where lines shunt parameters are neglecteq]] The method 
presented in fS\ always requires a base-case load flow solution 
and it relies on the assumption that all generators are PV nodes 
(i.e. with fixed voltage magnitude). Also, it does not account 
for the mutual coupling between different phase conductors. 
The approach that appears the more general among the 
above listed is the one proposed in lfT4ll . However, this method 
depends on a pseudo-load flow approach (i.e. it makes use 
of a Gauss-Seidel iterative process with a fixed number of 
iterations) which influences the accuracy of the computed 
coefficients. Furthermore, compared to lfT4l we have been 
able to; 

• generaUze the problem formulation for a generic number 
of slack busses; 

• extend the computation of sensitivities to tap-changers 
positions (i.e. changes of slack busses reference volt- 
ages); 

• provide the proof that the analytical computation of 
sensitivities admits a unique solution for the case of radial 
networks and 

• take into account the inherent multiphase and unbalanced 
nature of distribution networks. 

The structure of the paper is the following: Section |ll] 
focuses on the problem formulation by describing, in detail, 
the analytical procedure at the base of the proposed method. It 
also includes a proof of uniqueness of the solution of the linear 
system used to calculate the sensitivity coefficients for the case 
of radial networks. The same section also provides a computa- 
tional cost analysis of the proposed method versus traditional 



approaches. Section III validates the proposed method using 
the IEEE 13 and 34 node test feeders. Section HV] shows an 
application example of sensitivity coefficients related to the 
optimal voltage control in unbalanced distribution networks 
taking advantage of the possibility of computing them for 
all the phases. Section W\ provides the final remarks about 
possible applications of the proposed method. 

II. Problem Formulation 

A. Classical Computation of Sensitivity Coefficients in Power 
Networks 

In this paragraph we make reference to a balanced network 
composed by K busses. 

Traditionally, there are three proposed ways to calculate 
the sensitivity coefficients of our interest. The first method 
consists of estimating them by a series of load flow calcula- 
tions each performed for a small variation of a single control 
variable (i.e. nodal power injections, Pi,Qi) [4jj^ 



(1) 



dm 


Am 




dm\ 


A|/.,| 




dPi 


APi 




dPi 


APi 




dm 


Am 


APi ,^,=0 


dm\ 

dQi 


AQi 




dQi 


AQi 


AP,,.^i=0 
AQ..i^,=0 



It is important to observe that line shunt parameters are non-negligible in 
case of networks characterized by the presence of coaxial cables. These types 
of components are typical in the context of urban distribution networks 

-In the rest of the paper complex numbers are denoted with a bar above 
(e.g. E) and complex conjugates with a bar below (e.g. E). 



where Ei is the direct sequence phase-to-ground voltage of 
node i and 7^ is the direct sequence current flow between 
nodes i and j {i,j £ {1 • • • A'}). 

The second method uses the Newton Raphson formulation 
of the load flow calculation to directly infer the voltage 
sensitivity coefficients as submatrices of the inverted Jacobian 
matrix (e.g. H, 111, HIl, lEI, US): 



J = 



dP dP 

d\E\ ~de 

dQ dQ 



d\E\ de 



(2) 



It is worth observing that such a method does not allow to 

compute the sensitivities against the transformers tap-changers 

dQ 
positions. Additionally, as known, the submatrix _ is 

usually adopted to express voltage variations as a function of 
reactive power injections when the ratio of longitudinal line 
resistance versus reactance is negligible. It is worth noting 
that such an assumption is no longer applicable to distribution 
systems that require in addition to take into account active 
power injections. 

A third method is derived from circuit theory. In this 
method Tellegen's theorem is applied in power networks and 
the computation of sensitivities relies on the concept of the 
so-called adjoint networks (e.g. O, lUFl, fTf], HI, lfT9l ). 
This approach requires a base-case load flow solution in order 
to build a specific adjoint network that needs to be solved in 
order to infer the desired sensitivities. 



B. Analytical Derivation of Voltage and Current Sensitivity 
Coefficients 

This subsection contains the main analytical development 
of this paper related to the derivation of the voltage sensitivity 
coefficients |j 

1) Voltage Sensitivity Coefficients: the analysis starts with 
the voltage sensitivity coefficients. To this end, we derive 
mathematical expressions that link bus voltages to bus active 
and reactive power injections. For this purpose, a A"-bus 3- 
phase generic electrical network is considered. The foUowing 
analysis treats each phase of the network separately and, thus, 
it can be applied to unbalanced networks. 

As known, the equations that link the voltage of each phase 
of the busses to the corresponding injected current are in total 
M = SAT and they are given by: 
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abc 
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[E. 



abc 



(3) 



where [labc] = [/iJiJ^.-Jf ,/f ,/f ]^, [Eabc] = 
[ElEl,El...,E^,E^,E^]J. We denoted by a, b, c the 
three network phases. The [Yabc] matrix is formed by using 
the so-called compound admittance matrix (e.g. ||20l ) as 



As shown in subsection 



II-B2 



wardly derived from the voltage ones 



the current sensitivities can be straightfor- 



follows: 



where it has been taken into account that: 
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In order to simplify the notation, in what follows we will 
assume the following correspondences: [labc] — [h, ■■■,Im]'^, 
[Khc] = [El,.. .,EMf and 



[YabcJ 



Yn 



Yi 



M 



Yi 



M 



Ymm 



For the rest of the analysis we will consider the network 
as composed by S slack busses and N busses with PQ 
injections, (i.e. {1, 2, • • • Af } = 5 U TV, with 5 n TV = 0). 
The PQ injections are considered constant and independent 
of the voltage. In this respect, we are assuming that for 
each separate perturbation of nodal power injections, the other 
loads/generators do not change their power set points. There- 
fore, the computation of the sensitivities inherently accounts 
for the whole response of the network in terms of variation 
of both active and reactive power flows. Such a consequence 
allows to compute the sensitivities in the close vicinity of the 
network state. 

The link between power injections and bus voltages reads: 



jesuAf 



(4) 



The derived system of equations Q holds for all the phases 
of each bus of the network. Since the objective is to calculate 
the partial derivatives of the voltage magnitude over the active 
and reactive power injected in the other busses, we have to 
consider separately the slack bus of the system. As known, 
the assumptions for the slack bus equations are to keep its 
voltage constant and equal to the network rated value, by also 
fixing its phase equal to zero. Hence, for the three phases of 
the slack bus, it holds that: 



dp 



= yieS. 



(5) 



At this point, by using equation Q as a starting point one 
can derive closed- form mathematical expressions to define 
and quantify voltage sensitivity coefficients with respect to 
active and reactive power variations in correspondence of 
the N busses of the network. To derive voltage sensitivity 
coefficients, the partial derivatives of the voltages with respect 
to the active and reactive power Pi and Qi of a bus I G TV 
have to be computed. The partial derivatives with respect to 
active power satisfy the following system of equations: 
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(6) 
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dp 
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The system of equations (J6]) is not linear over complex 
numbers, but it is linear with respect to f^,gp^, therefore 
it is linear over real numbers with respect to rectangular 
coordinates. As we show next, it has a unique solution and 
can therefore be used to compute the partial derivatives in 
rectangular coordinates to reduce the computational effort. 

A similar system of equations holds for the sensitivity 
coefficients with respect to the injected reactive power Qi. 
With the same reasoning, by taking into account that: 





OQi 


dQi 


obtain that: 
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(8) 



(9) 



By observing the above linear systems of equations (|6| and 
(|9|, we can see that the matrix that needs to be inverted in 
order to solve the system is fixed independently of the power 
of the ^th bus with respect to which we want to compute the 
partial derivatives. The only element that changes is the left 
hand side of the equations. 

Once §p-,75p^ are obtained, the partial derivatives of the 
voltage magnitude can be expressed as: 



dm 

dp 
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Re{E, 



dEk 
dP, 



(10) 



and similar equations hold for derivatives with respect to 
reactive power injections. 

Theorem 1: The system of equations (J6]), where I is fixed 
and the unknowns are ^p^, i e TV, has a unique solution for 
every radial electrical network. The same holds for the system 
of equations ¥)\, where the unknowns are |^, i e TV. 

Proof: Smce the system is linear with respect to rectangu- 
lar coordinates and there are as many unknowns as equations, 
the theorem is equivalent to showing that the corresponding 
homogeneous system of equations has only the trivial solution. 
The homogeneous system can be written as: 



= A, 



jGSUAA 



Y,,E, 



E, E ^^I'^i ' 



yieAf (11) 



where A, are the unknown complex numbers, defined for i € 
TV. We want to show that A^ = for all i e TV. Let us 
consider two electrical networks with the same topology, i.e. 
same [Yabc] matrix, where the voltages are given. In the first 
network, the voltages are 



(12) 



El — E, , 


\/ieS 


El = E, + \ , 


\/i£N 


nd network they are 
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E'' = ^, - A, 
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(13) 



Let S'^ be the conjugate of the absorbed/injected power at the 
ith bus in the first network, and S" in the second. Apply 
equation (HI) to bus i E Af in the first network; 

s[ = E', E %,^ 

jesuAf 

= (E, + A,) ( E %,E, + E Y,,{E, + A,) I 
\je_s jeAf _ J 

^ E^ E Y^jE, + A, E Y,,A, 

+ A, E Y,jE,+E^ E %,A, 
jesuAf jeM 

Similarly, for the second network and for all busses i £ Af: 

S'l = E^ E Y^E, + A, E Y,,A, 

- A, E Y,,E, - E, E y„-A, 
jesuN' j€Af 

Subtract the last two equations and obtain 



s'-s'! = 



A. E 

jesuAf 



Y E- 



E, E Y,,K, 

jeAf 



By equation ( 1 1 1, it follows that S'l = S'l for all i E J\f. Thus 



the two networks have the same active and reactive powers at 
all non slack busses and the same voltages at all slack busses. 
As discussed in |21| for radial distribution networks such an 
assumption means that the load flow problem always has a 
unique solution. Therefore, it follows that the voltage profile 
of these networks must be exactly the same, i.e. _E, — A^ — 
Ei + Ai for alH e TV and thus A^ = for aU i E Af. ■ 

2) Current Sensitivity Coefficients: From the previous ana- 
lysis, the sensitivity coefficients linking the power injections 
to the voltage variations are known. Thus, it is straightforward 
to express the branch current sensitivities with respect to the 
same power injections. Assuming to represent the lines that 
compose the network by means of tt models, the current flow 
lij between nodes i and j can be expressed as a function of the 
phase-to-ground voltages of the relevant i, j nodes as follows: 



h,=%,{E,~E,) 



(14) 



where Yij is the generic element of [Yabc] matrix between 
node i and node j. 

Since the voltages can be expressed as a function of the 
power injections into the network busses, the partial derivatives 
of the current with respect to the active and reactive power 
injections in the network can be expressed as; 



dPi 
dljj 
dQi 
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(15) 



Applying the same reasoning as earlier, the branch current 
sensitivity coefficients with respect to an active power Pi can 
be computed using the following expressions: 



d\U 



dp 



i^.r^'^-^^ dPi 



(16) 



Similar expressions can be derived for the current coefficients 
with respect to the reactive power in the busses as: 



d\U 



dQi 






(17) 



C. Sensitivity Coefficients with respect to tap positions of 
transformers 

This subsection is devoted to the derivation of analytical 
expressions for the voltage sensitivity coefficientaj with re- 
spect to tap positions of a transformer We assume that 
transformers tap-changers are located in correspondence of 
the slack busses of the network as for distribution networks 
these represent the connections to external transmission or 
sub-transmission networks. As a consequence, the voltage 
sensitivities as a function of the tap positions are equivalent 
to the voltage sensitivities as a function of the slack reference 
voltage. We assume that the transformers voltage variations 
due to tap position changes are small enough so that the 
partial derivatives considered in the following analysis are 
meaningful. Furthermore, we assume that the power injections 
at the network busses are constant and independent of the 
voltage. 



With the same reasoning as in Sec. II-B[ the analysis starts 



in equation Q. We write E( = \Ee\e^^for all busses i. For 
a bus i E Af the partial derivatives with respect to the voltage 
magnitude \E]k\ of a slack bus fc G 5 are considered: 
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(18) 



where 
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We have taken into account that: 
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(19) 



(20) 



The derived system of equations (18i is linear with respect 
to Wik ^iid Wik, and has the same associated matrix as the 
system in (|6]l. Since the resulting homogeneous system of 
equations is identical to the one in ([TTJ, by Theorem [T] it has 
a unique solution. 

After resolution of ( [T8] l, we find that the sensitivity coef- 
ficients with respect to the tap position of the transformer at 
bus k are given by 



d\Ek 
d\Ek 



= \E,\Re 



E, 



(21) 



II-B2 



^Note as shown in Sec. 
the ones of currents can be computed directly. 



once the voltage sensitivities are obtained 



D. Computational Cost Analysis for Voltage Sensitivities with 
respect to PQ injections 

The aim of this subsection is to show the computational 
advantage of the proposed method compared to the classical 
approach with respect to the computation of voltage sensiti- 
vities as a function of power injections onljrl Furthermore, 
the two methods are applied to the IEEE 13 and 34 node 
test feeders and compared in terms of CPU time necessary to 
calculate the voltage sensitivity coefficients. 

We are assuming that: 

1) there are loads/injections in all three phases of the 
system and 

2) the phasors of phase-to-ground voltages in all the net- 
work are known (e.g. coming from a state estimation 
process {2T\). 

In the following table. Algorithm 1 shows the steps required 
to calculate the voltage sensitivity coefficients using the tradi- 
tional method and Algorithm 2 shows the corresponding steps 
using the analytical method proposed here. 

For the traditional method an updated Jacobian needs to 
be built, and its inverse will provide the desired voltage 
sensitivities. For the analytical method the corresponding 
steps refer to invert a square matrix of size 2N (as reported 
in Section III-Bll N refers to the number of network busses 
with PQ injections) and multiply the inverse matrix with one 
column vector for each PQ bus in the network. 

Algorithm 1 Computation of voltage sensitivity coefficients 
using the Jacobian method 
1: build Jacobian matrix associated to the Newton Raphson 

method 
2: invert matrix J of size 2A^ x 2A^ 

3: extract the sub-matrices corresponding to the desired sen- 
sitivity coefficients 

Algorithm 2 Computation of voltage sensitivity coefficients 
using the analytical method 

1: build the matrix of the linear system of equations 

2: invert matrix of size 27V x 2N 

3: do A^ multiplications of the inverse matrix with vectors of 
size 2A^ X 1 

In Table 11] the mean CPU time necessary to calculate the 
voltage sensitivity coefficients is presented for the IEEE 13 
and 34 node test feeders respectively, when 1000 iterations of 
the method are executed. It can be observed that the analytical 
approach exhibits an improvement of performance which is 
of 2.34 for the IEEE 13 node test feeder and 2.52 for the 
IEEE 34 node test feeder. In the same table the relevant 
95% confidence intervals are also reported for the computation 
of the coefficients for the two benchmark feeders. One can 
observe the advantage of the proposed analytical method as 
the number of busses in the network increases. It is worth 
observing that such an improvement depends not only on 



As already pointed out in Sec II-A traditional Jacobian based sensitivity 



the number of busses but also on the network topology (i.e. 
sparsity of the [Y] admittance matrix). 

Table I 

CPU TIME NECESSARY FOR CALCULATING VOLTAGE SENSITIVITY 

COEFFICIENTS IN THE IEEE 13 AND THE 34 NODE TEST FEEDERS WHEN 

ALL PHASES OF ALL BUSSES HAVE LOADS 



13 bus feeder 



Jacobian 

28.8 ±0.18 msec 



Analytical 

12.5 ± 0.43 msec 



ratio 

2.34 



34 bus feeder 209.8 ± 1.30 msec 83.4 ± 0.59 msec 2.52 



III. Numerical validation 

The numerical validation of the proposed method for the 
computation of voltage/current sensitivities is performed with 
two different approaches. In particular, as the inverse of the 
load flow Jacobian matrix provides the voltage sensitivities, 
the comparison reported below makes reference to such a 
method for the voltage sensitivities only. On the contrary, as 
the inverse of the load flow Jacobian matrix does not provide 
current sensitivity coefficients, their accuracy is evaluated by 
using a numerical approach where the load flow problem is 
solved by applying small injection perturbations into a given 
network (see Section |II-A[ ). A similar approach is deployed 
to validate the sensitivities with respect to tap positions of the 
transformers, i.e. small perturbations of the voltage magnitude 
of one phase of the slack bus and solution of the load 
flow problem. FiglT] shows the IEEE 13 nodes test feeder 
implemented in the EMTP-RV simulation environment (f23l. 
Ell, lEl) adopted to perform the multiphase load flow. 
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Figure 1. IEEE 13 node test feeder represented in the EMTP-RV simulation 
environment. 
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For the sake of brevity we limit the validation of the 
proposed method to a reduced number of busses exhibiting 
the largest voltage sensitivity against PQ load/injections. In 
particular, we refer to the variation of voltages at bus 8 with 
respect to load/injection in bus 9, i.e. 



d\EI\ d\El\ d\Ei\ d\El 



dp^ dp^ dQl dQi 



In Fig |2(a)| the voltage sensitivity of phase b bus 8 is shown 
with respect to active power absorption and generation at phase 
b of bus 9. We assume the convention that positive values of 
P and Q denote power absorption, whereas negative values 



correspond to power generation. Fig 2(b) shows for the same 
busses as Fig |2(a)| the same sensitivity but referring to voltage 
and power belonging to different phases. Additionally, Fig. 
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2(c) and 2(d) show the voltage sensitivity of bus 8 with respect 
to reactive power absorption and generation at bus 9. In all 
these four figures the dashed line represents the relative error 
between the traditional approach (i.e. based on the inverse 
of the Jacobian matrix) and the analytical method proposed 
here. As it can be observed, the overall errors are in the order 
of magnitude of 10^*'. In Fig 3(a) and Fig 3(b) the current 
sensitivity coefficient of phase a of branch 10 — 13 is presented 
with respect to active and reactive power absorption/generation 
at phase a of bus 13. In the same figures, the dashed lines 
represent the relative error between the analytical values and 
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Figure 2. Voltage sensitivity coefficient of phase a and h of bus 8 with 
respect to active and reactive power generation/absorption at phase b of bus 
9. 



Figure 3. Cun'ent sensitivity coefficients of phase a of branch 10-13 with 
respect to power generation/absorption at phase a of node 13. 



the numerical ones. Even for these coefficients extremely low 
errors are obtained. 

Concerning the validation of voltage sensitivities against 
tap-changer positions, we have made reference to the IEEE 13 
node test feeder where the slack bus and therefore the primary 
substation transformer is placed in correspondence of node 1 . 
We assume to vary the slack bus voltage of ±6% over 72 
tap positions (where position "0" refers to the network rated 
voltage). In Fig. l4] the sensitivity of voltage in phase a of 
bus 7 is shown w.rt. the tap positions in phase a, b and c 
of the slack. Also, in this case the difference between the 
analytically inferred sensitivities and the numerical computed 
ones is negligible (i.e. in the order of magnitude of 10^"*). 
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(b) Voltage sensitivity coefficient of phase a of bus 7 with respect to 
transformer's tap position at phase b of the slack bus. 
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(c) Voltage sensitivity coefficients 

to reactive power absorption at phase a of node 13 as 

a function of the distance from the slack bus. 
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Figure 4. Voltage sensitivity coefficient of phase a of bus 7 with respect to 
transformer's tap positions. 
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Figure 5. Voltage sensitivity coefficients with respect to power absorption 
at phase a of bus 13 as a function of the distance from the slack bus. 



It is worth observing that for the case of the voltage 
sensitivities, coefficients that refer to the voltage variation 
as a function of a perturbation (power injection or tap- 
changer position) of the same phase, show the largest coupling 
although a non-negligible cross dependency can be observed 
between different phases. 

Finally, FiglS] depicts the variation of voltage sensitivity 
coefficients in all the network with respect to active and 
reactive power absorption at phase a of bus 13 as a function 
of the distance from the slack bus in feet. 

This type of representation allows to observe the over- 
all network behavior against specific PQ busses absorp- 
tions/injections. In particular, we can see that larger sensiti- 
vities are observed when the distance between the considered 
voltage and the slack bus increases. Furthermore, a lower, but 
quantified dependency between coefficients related to different 
phases, can be observed. Also, as expected, reactive power 
has a larger influence on voltage variations although the active 
power exhibits a non negligible influence. 

From the operational point of view it is worth observing 
that, figures as Fig|5] provide to network operators an im- 
mediate view of the response of the electrical network against 
specific loads/injections that could also be used for closed loop 
control or contingency analysis. 
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Figure 6. IEEE 34 node test feeder represented in the EMTP-RV simulation 
environment. 



IV. Application of the proposed procedure to the 

PROBLEM OF OPTIMAL VOLTAGE CONTROL 

For the application part, the IEEE 34 test node feeder is 
considered as depicted in Fig|6] In busses 18, 23, 24 and 
33 we assume to have distributed energy resources that the 
Distribution Network Operator (DNO) can control in terms 
of active and reactive power. Their initial operating values, 
as well as their rated power outputs, are shown in Table [III 
Furthermore, the DNO has control on the transformer's tap 
positions. 

Table II 

Initial and maximum operational set points of the DERs and 

the tap-changers in the 34 test node feeder 





Pinit(kW) 


Pmax(kW) 


"init 


J^min 


Umax 


DERis 


210 


300 





-36 


+36 


DER23 


100 


600 


DER24 


250 


600 


DER33 


150 


300 



The optimal control problem is formulated as a linear 
one taking advantage of the voltage sensitivity coefficients. 
The controlled variables are the bus node voltages and the 
control variables are the active and reactive power injections 
of the DER and the transformer's tap positions under the 
control of the DNO, Ax = [APd^ja, AQ^i^;^, An]. It 
is important to state that, formally, this problem is a mixed 
integer optimization problem due to the tap positions of the 
transformers. However, for reasons of simplicity, the tap 
positions are considered pseudo-continuous variables which 
are rounded to the nearest integer once the optimal solution 
is reached. The objective of the linear optimization problem 
relevant to the problem is: 



mm 

Ax 



E, - E 



(22) 



The linearized relationship that links bus voltages with control 
variables is expressed in the following way (e.g. |4|): 



A|^,| = Kp.AP, + Kq^AQ, + K„,An, 



(23) 



where Kp, is the vector of sensitivity coefficients with respect 
to the active powers of the DERs, Kg is the vector of 
sensitivity coefficients with respect to the reactive powers of 
the DERs and K„j is the vector of sensitivity coefficients 
with respect to the transformer's tap positions. The imposed 
constraints on the operational points of the DERs and the tap 
positions are the following: 

< Pder, < Pder,,^^^ (24) 

QoERi^.^ < QoERi < QoERi^^^ 
'^min _ ^ _ ^max 

In order to simplify the analysis, we have assumed that the 
DER capability curves are rectangular ones in the PQ plane. 
The formulated linearized problem is solved by using the 
linear least squares method. The method used to calculate 
analytically the sensitivity coefficients allows us to consider 
two different optimization scenarios. In the first (opti), the 
operator of the system is assumed to control the set points 
of the DERs considering that they are injecting equal powers 



10 



into the three phases, whereas in the second case (opt2) it is 
assumed to have a more sophisticated control on each of the 
phases independently except for the tap-changers positions. 
It is worth noting that this second option, although far from 
a realistic implementation, allows us to show the capability 
of the proposed method to deal with the inherent unbalanced 
nature of distribution networks. Table |in] and Table |IV] show 
the optimal operational set points corresponding to these cases. 

Table III 

Optimal operational set points of the DERs and the 

tap-changers in the 34 test node feeder when the system 

operator has control on their 3-phase output 





Popti(kW) 


Qopti(kVar) 


Hopti 


DERis 


300 


300 


-2 


DER23 


600 


600 


DER24 


600 


264.06 


DER33 


300 


-14.46 



Additionally, in Fig|7] the voltage profile of the busses of 
the system is presented in the initial and the optimal cases. 
The solid line in the figures shows the initial voltage profile, 
the solid line with the markers shows the first case optimal 
scenario (opti) and the dashed line represents the second case 
where the DNO has full control in each of the phases of the 
DERS {opt2)- The offset in the graphs, observed in the slack 
bus, depicts the optimal tap position in each case. What can 
be observed is that, when there is a possibility to control each 
of the three phases of the DERs output, the optimal voltage 
profile is better than the one corresponding to control of the 
3 -phase output of the set points of the DERs. 

Table IV 

Optimal operational set points of the DERs and the 

tap-changers in the 34 test node feeder when the system 

operator has control on each of the three phases 

independently 





PoPt2(kW) 


Qopt2(kVar) 


"opta 


DERfg 


100 


100 


-Hi 


DERfg 


100 


-88.56 


DER-, 





83 


DERI3 


200 


200 


DER^a 


200 


200 


DER§3 





200 


DERI4 


200 


102.81 


DERg4 


196.51 


200 


DER-4 


111.40 


200 


DER« 


100 


-27.88 


DERga 


100 


100 


DER§3 


98.40 


100 



V. Conclusion 

In this paper we have proposed a new method for the 
analytical computation of voltages and currents sensitivity 
coefficients as a function of the nodal power injections. The 
contributions of the proposed method are the following: (i) it 
is generalized to account for a generic number of slack busses; 
(ii) it allows the computation of sensitivities w.rt. tap-changer 
positions (iii) it is proved to admit a unique solution for the 
case of radial networks and (iv) it supports the computation 
of the sensitivities for a generic unbalanced electrical network 
by using the [Y] compound matrix being, thus, suitable for 
distribution systems. 
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(a) IEEE 34 node test feeder - Voltage profile of phase a of the busses. 



1.02 

1.015 

1.01 

i' 1.005 

L 

: 1 
° 0.995 
i 0.99 

0.985 
0.98 

0.975 
0.97 



Phase b 




— initial voltage profile 

^^three phase balanced injections 

— three phase unbalanced injections 



50000 100000 150000 

Distance from slack bus (ft) 



200000 



(b) IEEE 34 node test feeder - Voltage profile of phase b of the busses. 
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(c) IEEE 34 node test feeder - Voltage profile of phase c of the busses. 



Figure 7. Initial and optimized voltage profile of the IEEE 34 node test 
feeder. 



Compared to the traditional use of the Jacobian load-flow 
matrix, it allows us to reduce the computation time by almost a 
factor of three, thus enabling, in principle, its implementation 
in real-time optimal controllers. 

The paper has also validated the proposed method by 
making reference to typical IEEE 13 and 34 nodes distribution 
test feeders. The former has been used to numerically validate 
the computation of the coefficients whilst the latter has been 
used to show an application example related to a possible 
integration of the proposed method for the problem of optimal 
voltage control in unbalanced distribution systems. 

It is worth observing that the proposed analytical computa- 
tion of voltages and currents sensitivities enables the reduction 
of the computational time of several traditional power systems 



problems involving non-negligible computational efforts, such 
as real-time centralized controls, contingency analysis or op- 
timal planning. 
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